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Abstract 

We prove a one-parameter family of sharp integral inequalities for 
functions on the n-dimensional unit ball. The inequalities are confor- 
mally invariant, and the sharp constants are attained for functions that 
are equivalent to a constant function under conformal transformations. 
As a limiting case, we obtain an inequality that generalizes Carleman's 
inequality for harmonic functions in the plane to poly-harmonic func- 
tions in higher dimensions. 



1 Introduction 

In this paper we compute the norm and the maximizers of some integral oper- 
ators between Lebesgue function spaces L p , used to solve the Dirichlet prob- 
lems for certain elliptic partial differential equations on the ball and the upper 
halfspace. One of the difficulties confronted is a lack of compactness due to 
conformal invariance of the equations in question. In some limiting cases, we 
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are able to differentiate these inequalities with respect to the parameter 1/p 
to obtain nonlinear weighted exponential-type inequalities governing solutions 
and subsolutions of higher order linear elliptic problems on the ball, which im- 
prove on the maximum principle and give higher dimensional generalizations 
of Carleman's isoperimetric inequality for subharmonic functions on the unit 
disc. 

In the following B n denotes the n-dimensional unit ball in Euclidean space, 
IMUp(n) is the LP norm of function u defined on domain Q, \B n \ is the vol- 
ume of B n and c(n,a,p) is some constant which depends on n,a and p. The 
parameter a satisfies 2 — n < a < 1. Before giving the main theorems, we 
will give an interesting corollary for the reason that it is easy to state and it 
is clearly a natural generalization of Carleman's inequality (see (8) below) in 
four dimension. 

Corollary 1. For any u satisfying A 2 u < and — |^ < 1, where 7 is the 
outer nomal of dB^, we have 



The sharp constant is assumed by the solution of A 2 u = in B n with boundary 
values — |^ = 1 and u = on dB A . 

For a function defined on R n_1 (thought of as the boundary of the upper 
half-space R+), we define a poly-harmonic extension as follows: for (X, x n ) G 



Here the choice of integrand guarantees independence of P a l on (X, x n ), while 
a < 1 ensures P a l < 00, and the normalization constants d r ^ a are chosen so 
that P a l = 1 (and can be expressed explicitly using T functions). Recalling 
that inversion in the unit sphere maps the halfspace x n > 1/2 to the unit ball 
centered at (0, 1), we see the conformal map 




(1) 



R" = R"" 1 x (0,+oo), 






(X, x n + %) 




(3) 
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maps the upper halfspace x n > to the standard ball : R" — > B n . 
Conformality of this map makes it easy to compute its Jacobian J(<f>) = 
\(X,x n + \)\~ 2n , and the Jacobian J(0| aR ») = \(X, |)|" 2(n_1) of its boundary 
trace. Indeed, pulls back the Euclidean metric g on B n to the conformally 
flat metric <p*g = \(X,x n + on R". Then it is not hard to check 

the formula 

f(X,x n ) = \(X,x n + l 1 )\ 2 - n - a fo ( j ) (X,x n ) (4) 

and its restriction to x n = boundary trace define Banach space isome- 
tries from / G L^+^(B n ) to / G L^V^(R™) and from L^=^ (dB n ) to 

2(n-l) , ~ ~ 

L"- 2 +«(R™ ) respectively. We define the poly-harmonic extension P a f of 

2(n-l) 

f g ^n-2+a (5_B n ) implicitly by using P a after pulling back from the ball to the 
halfspace: 

(PJ) o 0(X, x B ) = | (X, x n + \) r a ~ 2 P a (J^r^^j ■ (5) 

When a = 0, P a / again becomes the usual harmonic extension to the ball. 
Another case of special interest is a = 2 — n, in which case the conformal 
factors are suppressed so that Pi-^X — 1, and the isometric Banach spaces 
are both of L°° type. When n = 2k the extended function turns out to be k 
harmonic on the 2k dimensional ball, i.e. A fc P 2 -2fe/ — 0. 

2(n-l) 

Theorem 1. For any f E L n - 2 + a {dB n ), n > 2 and n — 2 + a > 0, we have 
the sharp inequality 

^h^(B^ s ^ f K*& t »„s (6) 

where the sharp constant S n ^ a depends only on n and a. The optimizers are 
unique up to a conformal transform and include the constant function f — 1 . 

We now study the limiting information. Letting f = 1 + "~ 2+a F and 
a — > 2 — n, we get the following inequality 

Theorem 2. For any F such that e F G L n ~ l (dB n ), n > 2, we have 

\\e In+p2 -" F \\LHB n ) < S n \\e F \\ L n- H9Bn): (7) 



3 



where I n = (log(X 2 + (x n + ±) 2 ) - d n;2 - n J Rn -i p^i^i log(F 2 + ^t^-^o 
C/p to a conformal transform any constant is an optimizer. 

In Theorem 2, when n = 4 we obtain Corollary 1. 

2 Some related literature 

Theorem 2 is a natural higher dimensional generalization of a well known 
inequality by Carleman [1] 

f e ^dx <-}-([ e u d6) 2 , (8) 

Jb 2 47T J dB2 

for all harmonic functions in B 2 . Equality occurs exactly for u = c and u = 
—2 log \x — Xo\ + c, where c is a constant and Xo G R 2 — B 2 . 

Although Carleman proved (8) initially for harmonic functions, it follows 
from the maximum principle that inequality (8) holds for subharmonic func- 
tions. Beckenbach and Rado [6] used Carleman's inequality to prove the 
isoperimetric inequality on a surface with non-positive Gauss curvature: If 
on a surface with non-positive Gauss curvature an analytic curve C of length 
L encloses a simply-connected domain D of area A, then the inequality 

L 2 > An A 

holds. This is exactly the sharp isoperimetric inequality in the plane. Their 
proof is quite simple: In isothermal coordinates (x, y) for a simply-connected 
domain D which is slightly larger than D, then the metric on D can be written 
as e 2w (dx 2 + dy 2 ), for (x,y) in some bounded domain f2 e R 2 . Now, the 
coordinate image of D in Q is a Jordan domain, so by the Riemann mapping 
theorem we can map it to B 2 conformally, which means D with the metric 
induced by the metric of the surface is isometric to (B 2 ,e 2u g), where u is a 
subharmonic function (By the non-positive curvature condition). Beckenbach 
and Rado's result now follows directly from Carleman's inequality. 

The 2008 paper by Hang, Wang and Yan [2], generalized this inequality to 
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higher dimensions as follows. For any harmonic function u in B n , 



\\U\\ _2n <Tl 2 (™- 1 )W„ 2n(n 11 2(n-l) , (9) 

where n > 3 and u n is the volume of B n . Any constant is an optimizer and 
it is unique up to a conformal transformation (as will be explained before 
the proof of Theorem 2). This is a special case of our Theorem 1 (a=0), and 
again because of the maximum principle, this inequality holds for subharmonic 
functions. Hang, Wang and Yan interpreted their inequality as the isoperimet- 

4 

ric inequality for B n with metric p"^g, where p is subharmonic (which means 
non-positive scalar curvature). By using the conformal map (2), the equivalent 
form of inequality (9) in the upper-half space is 



||Po/|| t5?W) <n~—oj n 2 ^ H/ll r ^i) w (10) 

2(ra-l) _, _, 

for all / G L ™- 2 (R n ). Here R™ is the boundary of R". The optimizers 

are f(Y) = , for some constant c, positive constant A and 

(A 2 +|Y-y | 2 )^- 

Y g R n_1 . 

Remark 1. We point out that the sharp inequality (9) combines with 
Brezis and Lieb's dual argument ([8] pagelO-11) to give the sharp version of 
inequality (1.9) in [8] when the domain is a ball: 

||V/|U2 (Bn) + C(n)\\f\\ 2^ > 5 B ||/|| - 

L «-2 (dB n ) L ° ^(B„) 

where S n is sharp Sobolev constant and C(n) can be determined by letting 
/ = 1 when the inequality becomes equality. This sharp Sobolev inequality 
with trace term was also proved by Maggi and Villani in [9] by using methods 
from optimal transportation. 

Remark 2. When — 1 < a < 1, from Caffarelli and Silvestre [4] we know 
u = P a f is the unique solution to the boundary value problem 

div(x a n Vu(X,x n )) = 0, (X, x n ) G R" 
U (X,0) = f,XeR n -\ 

Then the fractional Laplacian can be defined by using an analogue of the 
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Dirichlet to Neumann map (—A) 1 2 a f = — lim x a n u y . So, our equivalent form 



of inequality (6) on R™, namely \\P a f\\ L _ty_ < S n , a \\f\\j^± pro- 

vides a sharp estimate for the L n ~ 2 + a norm of solution of the above boundary 
value problem. 



3 Proof of Theorem 1 

Since P a enjoys very similar properties to the special case P (classical har- 
monic extension), we are also able to use the method of symmetrization de- 
veloped by Lieb [7] to prove the existence of maximizer as Hang, Wang and 
Yan did in [2]. The following Lemmas are parallel to those in [2], but notice 
that now we are dealing with poly-harmonic extension instead of harmonic 
extension. 

Recall if X is a measure space, p > and u is a measurable function on Q, 
then 

I 

||w|| L p = supt||w| > t\p. 

t>0 

The space L P W (X) = {u:u is measurable and < oo}. More generally, 

for any < p < oo and < q < oo, we have Lorentz norm || ■ \\ LP , g which 
is defined by ||w||li>,« = P*(f£°t q \\u\ > t|py)« and Lorentz space L p ' q (Q) . 
L^(f2) = L p '°°(£l) is a special case of such spaces. 

Lemma 1. Defining P a as in (1), there exist constants c n>a and c n>aiP such 
that 

\\Paf\\ ^ < C n , a ||/ 1| L i (Jr .-i) 

and 

\\ P *f\\ L ^ {R n +) < °n,a,p || /|| ^(B"- 1 ) 

for all 1 < p < oo. Moreover for 1 < p < oo we have 

Proof of Lemma 1. To prove the weak estimate, we may assume / > and 
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L i fR n-ix = 1. It is easy to see (P a f)(X,x n ) < -%^r for (X,x n ) G R+ and 
[ (P a f)(X,x n )d n ^Xdx n 

J(X,x n )£-R.l,0<y<b 

= [ d^Y{f{Y) fdy [ d n J;f — d^X) 

Jr.™- 1 Jo Jr" 1 ((X-Y) 2 + xl) — 



for b > 0. Hence for t > 0, 



|P„/>f| = |{(X,x n )GR^:0< a ; n <K»"^ T ,( J Pa/)(X, a ; n )>t}| 

< 7 / x {P a f){X,x n )d n - l Xdx n 

1 JxeR"- 1 ,o<i„<(t 1 a () 

The weak type inequality follows. The strong estimate follows from Marcinkiewicz 
interpolation theorem (see [5], pl97) and the basic fact ||-P a /|U°°(R™.) < II/IIl^r™- 1 )- 

Remark 3. In fact when p = 2 ^~^ and a = 0, the second estimate was 
also proved by Brezis and Lieb [8] by using some elementary dual argument. 

Lemma 2. If n > 2 and 1 < p < oo, then the supremum 

np np 

Cn7a] P = SU P {||P a /||"- : H/ll^fln-i) = 1}, (11) 

is attained by some function. After multiplying by a nonzero constant, every 
maximizer f is nonnegative, radially symmetric with respect to some point, 
strictly decreasing in the radial direction and it satisfies the following Euler- 
Lagrange equation 

forr x = I „ y SC 2 ^ {Paf)^-\x, Xn )d n - i xdx n . 

Jri {{X -Yf + xl) 2 

In particular, if n > 2, p = and n — 2 + a > 0, then every maximizer is 

of the form 

f(Y) = ±c(n, a)( x2+ Y_y (12) 
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for some X > 0, Y G R 



Proof of Lemma 2. First we recall the important Riesz rearrangement in- 
equality. Let u be a measurable function on R n , the symmetric rearrangement 
of u is the nonnegative lower semi-continuous radial decreasing function u* 
that has the same distribution as u. We have 

/ dx I u(x)v(y — x)w(y)dy < / dx u*(x)v*(y — x)w*(y)dy. 
Jr" Jn n Jn n Jr» 

Using the fact ||w||,Lp(R n ) = ||w*||lp(r") for p > and the standard duality 
argument, we see for 1 < p < oo, 

\\u * ^||lp(R") < \\u* * V*\\ L P(fi n )- 

Moreover if u is nonnegative radially symmetric and strictly decreasing in the 
radial direction, v is nonnegative, 1 < p < oo and 

||w * w||lp(R") = \\ u * * ' u *||lp(R' 1 ) < OO, 

then for some x G R™, we have v(x) = v*(x — x ). 

Now, assume f\ is a maximizing sequence in (11). Since ||/*||lp(r,™-1) = 
11/ilL^R"- 1 ) = 1 an d 



np I np 

I n-1 / 1 1 73 > -P \\ n ~ 1 



np — / ||-Pct,^ n * fi\\ np 1 dlj 

"l^(r^) L "l^cr"- 1 ) 



POO 

Jo 



np 

^ / -Pa,^n * fill -HP , ,dx n 



np 

IIP || 1 

L^T(R") 

where P a ,x n = ^n,a (xa n +a .s ) an d notice that it is symmetric and strictly de- 
creasing in the radial direction of X variable for any fixed x n . We see f\ is 
again a maximizing sequence. Hence we may assume fi is a nonnegative radial 
decreasing function. 

For any / G L p (K n ' 1 ) and any A > 0, we let f x (Y) = A~^/©, so that 
is clear that (P a f x )(X,x n ) = X^ (P a f){f, ?f) and hence H/i^n-i) = 
LP(R n-i) and HPo/*!!^^ = WPafW^^y For convenience, denote 



ei = (1,0, - •• ,0) e R and 

a, = sup{^( ei )|A > 0} = S up{A^/ 4 (^)|A > 0}. 

n-l 1. 

It follows that < < aj|y| " , and hence ||/i||.z>°°(R™-i) < |-B n _ 1 |pa i 

Now 



H P ^Hl^(r;) - c ( n ' a ^)II^H L -^ (R n-i ) 

n-l j_ 

< c(n,a,p)||/ i || I ^ (R „_ 1) ||/ i ||£ I)i00(R „_ 1) 



< c(n,a,p)a?, 



which implies a, L > c(n,a,p) > 0. We may choose Aj > such that f i *(ei) > 
c(n,a,p) > 0. Replacing /j by /j Al we may assume /(ei) > c(n,a,p) > 0. On 
the other hand , since fi is nonnegative radial decreasing and H/iHl^r™- 1 ) = 1) 
we see 

\fi(Y)\ < \B n ^\- l v\Y\-^ . 

Hence after passing to a subsequence, we may find a nonnegative radial de- 
creasing function / such that /j — >■ / a.e. It follows that f(Y) > c(n, a,p) > 
for \Y\ < 1, ft -± f in L^R™" 1 ) and ||/|| LP(R n-i) < 1. By Lieb [7](Lemma 
2.6), we have 

I \\h\ p -\f\ p -\h-f\ p \d n ' 1 y^^ 

It follows that 

ll/i — f\\'LP('R n - 1 ) = II^IIlp(R 11 - 1 ) ~~ II ^ II LP(R n_1 ) + °W 

= l-H/llip(R-l)+0(l)- 

On the other hand, since (P a fi)(X,x n ) — > (P a f)(X,x n )) for (X,x n ) e R™ and 
ll p a/i|| L ^ r(R „ ) < c w , we see 

np np np 

II^IR , = +11^-^/11%, +0(1) 

L^(RJ) L^(R^) (RIJ.) 

np np np np 

< c ",a!p 11/ 1| 2p(r™-1) + C n,a!p||/i ~~ / II + 
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np np 



rip 



Hence 

1< ||/||^ R „- 1) + ||/,-/||^ r R „- 1) + o(l). 
Let % — > oo, we see 

i^ii/ii^-ij + a-n/iiipfR-- 1 )^- 

Since^- > 1 and / ^ 0, we see ||/||lp(r«-1) = 1. Hence /; ->■ / in L P (R" _1 ) 
and / is a maximizer. This implies the existence of an extremal function. 

Assume / G ^(R™" 1 ) is a maximizer, then so is |/|. Hence ||-Po/|| L? ^ r , R „> = 
\\Pa\f\\\ L ^ {nn+ y On the other hand , since \(PJ)(X,x n )\ < (P a \f\)(X,x n ) 
for (X,x n ) G R™, we see \P a f\ = P a (\f\) an d this implies either / > or 
/ < 0. Assume / > 0, then the Euler-Lagrange equation after scaling by a 
positive constant is given by 

f{Yf^= I f~ a (P a f)^-\X,x n )<T-'Xdx n . 

Jul {{x - yy + x 2 n ) 2 

On the other hand, we know for x n > 0, 



\P ax *f\\jm r , = llPax * f ' 



np 



which implies f(Y) = f*(Y — Y ) for some Y . It follows from the above 
Euler-Lagrange equation and Lemma 2.2 of Lieb [7] that / must be strictly 
decreasing along the radial direction. 

2( n \) 2(n— 1) -J 

For the case when p = ^_ 2+a , we first observe that if / G L«-2+a (R n ~ ), let 
u = P a f, f = jyparz /(j^jg) and u = |( X;:rn ) 1 |n-2+a /( | ( ( ^ n n ) |2 ), then we have u = 

Paf, \\f\\ 2(n-l) = llfll 2(n-l) and llull 2n = llull 2n 

" J "i^fR"-!) L^TSCR™- 1 ) "l^T^+^(R") 11 M LS=2+S(R™) 

This is the conformal invariance property for the particular power. As a conse- 
quence, if / is a maximizer which is nonnegative and radial, then | y | J. 2+a f(j£^~ 
ei) is also a maximizer. In particular, f (j^s ~ e i) is radial with respect 

to some points. To find such /, we need the following useful Proposition of 
Hang, Wang and Yan [2]( Proposition 4.1). 

Lemma 3. Let n > 2, u be a function on R n which is radial with respect 
to the origin, < u(x) < oo for x ^ 0, e± — (1, 0, • • • , 0), a G R, a ^ 0. 
If v(x) = \x\ a u(j^2 — e i) is radial with respect to some point, then either 
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u(x) = (ci |x| 2 + C2) 2 for some C\ > 0, C2 > or 



u(x) = < \ T (13) 

C2 if x — 0. 



Proof of Lemma 2 continued. Since ||/|| 2(n-i) = 1 and it is strictly 

decreasing along the radial direction, we have < f(Y) < 00 for all Y ^ 0. 
Note that since / satisfies the Euler-Lagrange equation, it is defined every- 
where instead of almost everywhere. It follows from Lemma 3 that f(Y) = 
(cx|F| 2 + c 2 y" 2 +a for some Ci,c 2 > (since / can not be constant function 
and the scalar multiple of |F| 2 ~™ is ruled out by the integrability) . Using the 
condition ||/|| gjn-i) = 1, it is easy to see CiC 2 = c n ^ a . Hence for some 

A > 0, 

f(Y) = c(n,a)( x2 + l *_ Yol2 )^. 

~ 2(n-l) ~ 

Proof of Theorem 1. For any / e L«- 2 + a (dB n ), let u = P a f, 
f = |(F,0) + (0,|)|™- 2 +«^° ' 

and 

By definition (4) we have u = P a f and by the discussion below (3) we have 
11/11 2(n-i) = 11/11 2(n-i) and IImII 2 ? = 1 1 it 1 1 gn Then, 

Theorem 1 follows easily from the above facts and Lemma 2. 



4 Proof of Theorem 2 

First we will discuss some conformal invariance properties of the operator P a . 
Let t be a conformal transform from B n to itself, r = t\qb„ is the induced 
conformal transform from dB n to itself, J is the Jacobian of r , J is the Jacobian 

2(ra-l) 

of r, £ = n — 2 + a. For / e L™- 2 + a (dB n ), when e 7^ 0, we have 

P a (J^f or) = J^ (P a f) o r. (14) 
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It is straightforward to check this property by using the definition of P a in (5). 
Now, for smooth function /, when e goes to it is obvious that 

P 2 -n(/°r) = (P 2 _ n /)of. (15) 

By letting f = 1 and taking derivative with respect to e at 0, we have 

—\s=« + P^-ni^^Tj log J) = ~dT U ° r + Yn log J (16) 

So the inequality in the Theorem 2 is invariant when F is replaced by F o r + 
log J. 

n— 1 & 

Proo/ of Theorem 2. Recalling p2_„l = 1, let / = 1 + eF, where F is some 
smooth function defined on dB n . By Theorem 1, we have the inequality 

\\P a {l + eF)\\ 2. <S n , a \\l + sF\\ 2(«-i) 
which means 

(/ (P a l)¥(i + £^)^)^<4 a (/ (l + eF^dO^. 

J B n P a 1 J dB n 

Note that when F = the above inequality becomes equality, then by the fol- 
lowing estimates we will see in this case the integrals in both sides will converge 

2 

to some finite numbers, which means the constant Sn,a will also converge. 

In order to take limit e — > 0, we need to apply the Dominated Convergence 
Theorem, we will bound the term P a l from below by a constant A and bound 
(P a l)~r from above by a constant B, both A and B are independent of e. Let 
us derive the lower bound of P a l first. From (1) and (4) we know 

1 C 1 



2 Jr"- 1 ((X - F) 2 + a;2)— (F 2 + 4)2 
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by letting U = in the integral, we have 



P a l)o0 = dn, a — ' ! £ rf C/ 

Jr"- 1 (u 2 + 1)— ((t/:r n + X) 2 + i)2 



i 2 + l) 2 t a ((^ n + x) 2 + |)- 

> r/„ n / ^-r:-;;- — — — ' "J ^." r/" "'f 



1,0 / 
J\u 



1 (|X| 2 +(a; n + ^U 

1 1 



^(1^ + ( 2 |X| 2 + 2:r 2 + i)i 



= A. 

Then, lets derive the upper bound of (P a l)^r, it is enough to prove that 
(P a l)~ is bounded from above by some constant B independent of e. As in 
the proof of lower bound, after the same change of variable we have 

r>A = n(Y <r \\ J / 

R -i ( |f / |2 + 1) ^(( f/a;n + X) 2 + i)f 



{Pal)^o<t> = g(X,x n )[d n , a — ■ „-„„_■ -^W-T^ dn ~ 1U 

\ Jr" 1 {\u\ 2 + 1)— {(Ux n + xy + i )2 

< g(X,x n )d n , a [ — — 1 _ a 1 ^ (T^U 

Jr- 1 (|f/| 2 + 1)— {{Ux n + X) 2 + \)~ 

< g (x,x n )^^d n0 [ n jd n - l U 



g(X,x n ) dn f n d nfl I 

«n,0 JR"- 1 



4' 

1 d" L Y 



n-2 



d n ,2-n 
d n ,0 

B, 



where g(X ) x n ) = {\X\ 2 + (x n + |) 2 ) 2 . For the first inequality we applied 
Jensen's inequality, since d na 1 n _ a is a probability density in R n_1 and 

n — 2 

g(t) = t~^~ is convex when t > 0. The last identity holds because 



x n 1 

'" ° ((x-y) 2 + x 2 Shy 2 + ^ 



d n ' l Y 



is the harmonic extension of function (Y 2 + j) 2 which is easy to verify that 
it is exactly ^=5-. 

(X^+{x n + \) 2 )^- 
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Now we can take limit e — > safely. By denoting 



1 n — z 1 e=0 

de 



= fax* + (x. + i» - < 2 -„ I _ > og P' 2 + 

we get ||e /n+2P2 - nF ||Ln(B n ) < S^He^HL"-^^)- After replacing 2F with F, the 
inequality in Theorem 2 is proved. Since constant functions are optimizers for 
the above inequality, conformal invariance of the inequality tells us that the 
functions 



F = C + log J 

n — 1 



are also optimizers. 



Remark 4. The uniqueness is lost when taking limit in the proof of 
Theorem 2. It would be interesting to find a suitable method to prove that 
the optimizer is unique up to conformal transforms. The main difficulty seems 
to the author is that the integral kernel is too complicated. 



5 Proof of Corollary 1 

Now we are in the situation where n = 4 and a = —2. By [12] we know that 

3 

(P- 2 f)(X,x n ) = rf 4 „ 2 J R3 ((x_ Y )2 n +x 2)n-i f(Y)dY is the bi-harmonic extension 
of the function f(Y) with boundary condition a(P " 2 /^ X ' a " ) L=o = 0. It is 
straightforward to check that under the conformal map the bi-harmonic 
property and the Neumann boundary condition are preserved in dimension 
four, we have that P-29 is a bi-harmonic extension of a function g defined on 
S* 3 to a function on _B 4 with boundary condition 9P ~^ 9 \ y=Q = 0. In view of 
Theorem 2, in order to prove Corollary 1 we only need to verify I 2 satisfies 
A 2 / 2 = 0, I 2 = 0| 53 and = 1. 

From the formula for 7 n , we have 
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By using the explicit formula of <f> one can get 



h = 2 log \ V -S\- 2C J } vl J log |e - S\d£, 

where 77 is a point in _B 4 , £ is a point on S 3 , S is the south pole of S 3 and 
C is the normalizing constant such that C f s3 ^Z^Le d£ — 1- Now from [12] 
we have the representation formula for bi-harmonic functions, namely for a 
smooth bi-harmonic function g on B 4 we have 

*>W,Wrf o *W,W<->* (l7) 

where D is a known constant. Although the function log |?7 — S\ is singular 
at the south pole S, If we apply the forthcoming approximation process we 
obtain 



h = 2\og\r)-S\-2C \ log g- 



J53 r 



since 



■ dlog Q^~ S ^ (£) = — 1. In the above equality — D J s3 is the bi- 

harmonic extension of constant function with boundary condition — §^(£) = 
1, so I 2 satisfies all three conditions mentioned above. 

Since log \rj — S\ is singular, we use approximation to justify the previous 
formula for I 2 . Take a sequence 

S t = (0,0, 0,-t)^S = (0,0, 0,-1), 

as £ — > 1 + . Then log \rj — S t \ is a smooth bi-harmonic function on _B 4 , so we 
have 

here we use y to denote the last coordinate of £. For fixed 77 G _B 4 , when t 
approximates 1 from the right, | log |£ — < | log |£ — S\\ for £ in a small 
neighborhood of S, since | log |£ — S\\ is integrable on S* 3 , by the Dominated 
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Convergence Theorem 

Js* \v-£r J \v-£\ 6 

as t — ?■ 1+. Similarly, when t is close to 1 from right hand side, we have 

, l+ty 1 10 

'l + t 2 + 2ty' ~ 2 |£-S| 2 ' 

Since \ + igigp is integrable on S* 3 , by the Dominated Convergence Theorem 
again we have 

D / 1 + ^ 1 I (1 - Iff 

as t — ?• 1+. Now by taking limit t — > 1+, it is clear that we have the represen- 
tation formula for log 1 77 — jS' | . Finally since the kernels in the representation 
formula (17) are positive, we conclude that the inequality in Corollary 1 is true 
for sub-biharmonic function u with boundary conditions — |^ = 1 and u — 
on dB 4 . 
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